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COORDINATE SYSTEM INFLUENCE ON THE REGULARIZED 
TRAJECTORY OPTIMIZATION PROBLEM 


By J. M. Lewallen, Manned Spacecraft Center, 

and 0. A. Schwausch, Lockheed Electronics Company 

SUMMARY 

This investigation studies the effect of using regu- 
larized variables to enhance the numerical integration 
process associated with the optimal trajectory of a con- 
tinuously thrusting space vehicle. The integration charac- 
teristics^of both the rectangular Cartesian and polar 
cylindrical'' coordinates are considered for an optimal, low- 
thrust, Earth-escape, spiral trajectory. The numerical 
accuracy achieved and the computer time required are compared 
for various numerical integration error bounds, by using 
both the unregularized and regularized equations. The results 
obtained indicate that for space vehicles which experience 
wide variations in the gravitational force magnitude, signi- 
ficant reductions in computing time can be obtained by 
using the regularized trajectory optimization equations. In 
some cases, the computing time is reduced by a factor of 
three if regularized variables are used. Furthermore, for 
the problem considered here, use of the polar coordinates 
consistently results in more favorable computer times than 
when rectangular coordinates are used. In addition, if the 
numerically evaluated Hamiltonian, which is theoretically 
constant, is used as an indication of integration error 
generation, the trade-off between integration time and inte- 
gration error becomes apparent. Finally, it is shown that 
the polar coordinates are less sensitive than the rectangular 
coordinates to errors in the initial Lagrange multipliers. 
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INTRODUCTION 


During the past decade, considerable effort has been 
directed toward determining numerical methods for optimiza- 
tion of nonlinear, dynamic systems. A comparison of the 
characteris tics of several of the more popular direct and 
indirect numerical optimization methods is given in Ref. 1. 
Further investigations dealing with the procedures for 
accelerating convergence of the indirect optimization 
methods are discussed in Ref. 2. The primary consideration 
in evaluating an optimization method is the computing time 
required for convergence to a sufficiently accurate solu- 
tion. These characteristics may be influenced by the func- 
tional form of the equations of motion as well as the choice 
of the coordinate system in which the motion is computed. 

Regularizing transformations have been used in celes- 
tial mechanics to eliminate singularities associated xvith 
gravitational force centers. Results reported in Ref. 3 
indicate that the numerical integration characteristics can 
be enhanced considerably when a regularized set of differen- 
tial equations are used for trajectories that experience 
close primary body approaches. This conclusion has been 
reached also in Ref. 4 for a wide range of problems in 
celestial mechanics. Based on these conclusions, a study 
was made of the applicability of using regularizing trans- 
formations to the problem of improving the computational 
characteristics of numerical optimization procedures. The 
results described in Ref. 5 indicate significant numerical 
advantages in terms of . computational time and accuracy of 
terminal condition satisfaction if regular variables are 
used. 
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The effect of the regularising transformation is 
obviously dependent on the choice of the coordinate system 
for the unregularized variables. The influence of the coor- 
dinate system on numerical error generation in the two-body 
problem has been studied in Ref. 6 and in the unregularized 
trajectory optimization problem in Refs. 7 and 8. These 
investigations indicate that the coordinate sytem used can 
have a significant effect on computation time and the accuracy 
of the resulting numerical solution. In particular, these 
investigations revealed that the polar coordinates were 
computationally superior to the rectangular coordinates for 
the continuously powered escape spiral. 

In the investigation discussed in the following section, 
the effect of using both rectangular Cartesian and polar 
cylindrical coordinate systems is studied for a minimum time, 
l,ow- thrust, Earth escape spiral. The numerical accuracy, 
the computation time and the convergence characteristics are 
compared by using both the regularized and unregularized 
equations for various bounds on the integration error. 

FORMULATION 


If the transfer trajectory for a continuously powered, 
low-thrust, space vehicle is to be time optimal, the following 
equations must be satisfied in the interval t t <_ t : 
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The quantity m = m 0 - gt where 3 is a constant mass flow 
rate and A and to are Lagrange multiplier vectors. The 
boundary conditions that must be satisfied are 
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By using a generalization of the classical Sundman regu- 
larizing transformation discussed in Ref, 9, i.e., 

, - 3/2 

dr = r dt (6) 

a set of regularized equations for the optimal trajectory 
can be obtained as follows: 
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where the primes indicate derivatives with respect to the 
pseudo time variable t rather than the real time t . 

This transformation is discussed in Ref. 5 where it is 
shown that Eqs . (7) and (8) are mathematically regular. This 
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vector form of the regularized equations is invariant with 
the choice of coordinate system. Hence, Eqs. (1) and (2) 
describe the optimal process in the unregularized rectan- 
gular and polar coordinates, while Eqs. (7) and (8) describe 
the regularized equations associated with each of the coor- 
dinate systems. Either set of equations represents the 
usual, nonlinear, two-point boundary value problem. 

DISCUSSION OF RESULTS 

From the preceding section, it is seen that the solution 
to the optimal trajectory problem involves the solution of a 
nonlinear two-point boundary value problem. Usually, efforts 
are made to obtain a numerical solution to Eqs. (1) and (2) 
which satisfies the boundary conditions given by Eqs. (3] , 
(4], and (5). Since Eqs. (3) specify only half the neces- 
sary initial conditions, values for the remaining unknown 
initial conditions, usually Lagrange multipliers and the 
unknown time, must be assumed before a numerical solution 
can be determined. Inasmuch as the values of the unknown 
initial boundary conditions are arbitrarily selected, the 
terminal constraints given by Eqs. (4) and (5) will not be 
satisfied. These arbitrarily selected initial conditions 
are changed systematically on subsequent iterations until 
the terminal constraints are satisfied more exactly. There 
are numerous procedures for obtaining the corrections to the 
unknown conditions. Several of the currently popular itera- 
tion procedures are discussed in Ref. 1. 

Adequate satisfaction of the specified terminal con- 
straints as well as sufficient numerical accuracy, must be 
achieved if an acceptable numerical solution is to be 
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obtained. Adequate terminal constraint satisfaction is 
obtained by requiring the norm of the terminal constraint 
error to be less than 10 . Sufficient numerical accuracy 

is obtained by using full-double precision arithmetic on 
the UNIVAC 1108 at the NASA Manned Spacecraft Center and 
by perform-ing the integrations with a variable step-size 
integration scheme, thereby maintaining the single-step error 
within certain desired tolerances. The integration scheme 
employed is a modified version of the scheme discussed in 
Ref. 10. This scheme uses a fourth-order Runge-Kutta 
starter and a fourth-order Adams -Bashford predictor corrector. 

In order to determine the individual effects of the 

coordinate system and regularization, the same problem must 

be solved in both coordinate systems and in both unregu- 

larized and regularized form. The optimal Earth escape 

spiral for a low-thrust space vehicle is an excellent 

example problem for regularization investigations since the 

2 

gravitational force magnitude varies by approximately 10 , 
and hence it is expected that a wide range of numerical 
integration step sizes will be required to maintain certain 
specified error bounds. 

Figure 1 shows the optimal escape spiral. Initially, 
the spacecraft is in a circular near-Earth orbit with a 
radius equal to 1.05 times the Earth radius. For a constant 
low- thrust space vehicle subjected to a thrust to mass ratio’ 
of 0.1, the spacecraft acquires escape energy in approxi- 
mately 70 normalized time units (approximately 15.7 hours) 
and reaches an orbit of radius equal to 8.5 times the Earth 
radius. Although this thrust to mass ratio is relatively 
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large, it was selected to compromise between a computationally 
expensive realistic trajectory and an inexpensive unrealistic 
one. The trend of the results is probably unaltered. Figure 
1 also shows the optimal control programs for both the rec- 
tangular and polar coordinate systems. Figure 2 shows the 
relationship between the real and regularized time for the 
optimal trajectory. 

Tables 1 through 3 compare the integration characteris- 
tics of the regularized and unregularized polar and rectan- 
gular coordinate systems for various absolute single-step 
integration error bounds. The error-bound separations in 
Tables 1, 2, and 3 are 10 6 , 10 4 , and 10 2 , respectively. 

The numerical integration characteristics which are compared 
include the amount of computer time needed to perform all 
integrations for the final converged iteration, the average 
amount of computer time required per integration step, the 
number of integration steps required, the number of step size 
changes made, and the norm of the terminal constraint error. 

The integration time shown in Tables 1 through 3 
represents the computation time needed to integrate the 
state equations, the Euler-Lagrange equations, and the 
perturbation equations from the initial time to the final 
time. The values shown also include the time required to 
monitor the single-step integration error and determine 
the appropriate integration step size. The appropriate step 
size is determined by comparing the single-step error with 
the desired accuracy limits. If either the maximum or 
minimum error limit is encountered, the step size is either 
halved or doubled. If, by doubling the step size, the maximum 
bound is violated, then the step size remains unchanged. The 
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total number of integration steps taken in the interval and 
the number of step-size changes necessary to maintain the 
desired accuracy are recorded also. No distinction is made 
in the Tables between step-size changes associated with 
doubling and halving. The average computer time per inte- 
gration step is recorded to indicate the degree of complexity 
of the equations for each case. Finally, in order to indicate 
the degree to which the terminal constraints are satisfied, 
the norm of the constraint error is recorded. This quantity 
should be considered with some reservation since the routine 
simply requires that the norm be less than 10 -7 . The extent 
to which this criterion is exceeded is not controlled and is 
an indication of the convergence rate. However, it also 
depends on how close the terminal norm for the previous 
iteration was to the required value of 10 -7 . 

The results presented in' Table 1 are for the relatively 
large error-bound separation of 10^. It is seen that the 
regularized variables, in either coordinate system, require 
considerably less computation time per iteration than the 
unregularized variables. In some cases, the time is reduced 
by a factor of three. The reason for the large saving in 
time is readily apparent when the combination of time per 
iteration step and the total number of steps is examined. 
Although the regularized equations are more time consuming 
to evaluate, as indicated by the time required per step, the 
large number of steps taken by the unregularized system of 
equations quickly causes the total time to exceed that of the 
regularized systems. Table 1 also indicates that the polar 
coordinates generally require less computer time than the 
rectangular coordinates. 
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The results shown in Table 2 for an error-bound separa- 
4 

tion of 10 agree with those presented in Table 1 and sub- 
stantiate the previous conclusions. Again, the regularized 
variables require less total computer time than the unregu- 
larized variables, and the polar coordinate systems exhibit 
shorter integration times than the rectangular coordinate 
systems. However, for this error-bound separation, the 
computation time advantage of the regularized systems has 
been reduced slightly. Note also that the difference in the 
total number of integration steps beti^een the regularized 
and unregularized variables has been reduced. In addition, 
the number of step-size changes for the regularized variables 
is less than the number of changes required by the unregu- 
larized variables. This is in keeping with the regulariza- 
tion theory which predicts that regularized variables will 
undergo fewer step-size changes than unregularized variables, 
provided a certain integration accuracy is to be maintained. 
(For the previous error-bound separation of 10 6 , a comparison 
of the number of step-size changes is invalid since, in some 
instances, the lower error bound was never encountered.) 

The results presented in Table 3 for the error-bound 
2 

separation of 10 generally agree with the results of Tables 
1 and 2. As in the previous tables, the polar coordinate 
system requires shorter integration times than the rectan- 
gular system. However, for this magnitude of error-bound 
separation, the integration times for the regularized and 
unregularized variables are essentially the same. The 
departures from the previously indicated trend can be 
explained by examining Table 4. 
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Shown in Table 4 are the error-bound encounters for 
certain integration error tolerances. The top line in each 
set of four lines represents the upper or maximum allowable 
error bound. Each succeeding line represents the minimum 
allowable error for a particular error-bound separation. 
Thus, the first set of four lines represents the integration 
error bounds of 10~ 4 and 10“ 6 , 10~ 4 and 10“ 8 , and 10“ 4 and 
10 . The boundary encounters are plotted as a function 

of the normalized trajectory time. One of the appropriate 
symbols, keyed in Table 4, records the encounter of the 
numerical error magnitude with either of the boundaries. 

An encounter' with the lower bound means the step size will 
be doubled; an encounter with the upper bound means the step 
size will be halved. 

Table 4 indicates that by maintaining the small inte- 
gration error-bound separation of 10 2 , the error in the 
unregularized rectangular variables is such that the step 
size is doubled three times during the escape trajectory 
for the 10 4 to 10 6 accuracy limits. Upon increasing the 
error separation to 10 4 , to give error bounds 10 -4 to 10” 8 , 
the unregularized rectangular error becomes less than the 
minimum acceptable error only twice, with the first boundary 
encounter coming after the 10 6 bound in the previous case 
had already been crossed twice. By doubling the step size 
early in the trajectory flight time in the 10 4 to 10~ 6 
case, 7 seconds of computer time were saved per iteration. 
This time saving was increased to approximately 10 seconds 
when comparing with the 10 -4 to 10 _1 ° accuracy level since 
the lower boundary for this case was never encountered. 

Thus, by requiring the rectangular error to be within the 
10 4 — 10 6 accuracy level rather than the 10 -4 — 10 -8 
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accuracy level, 253 integration steps were eliminated. 
Elimination of these 253 steps, each consuming approximately 
.0276 seconds of computer time, resulted in saving 7 seconds 
of computer time per iteration. Likewise, by requiring the 
integration error to be within the 10 4 — 10 accuracy level 
rather than the 10 4 — 10 10 interval, a 10 -second saving 
in computer time per iteration was realized. This same trend 
appeared in both the rectangular and polar coordinates, for 
the other error bounds shown. By maintaining the integration 
error within the smaller error bounds, the total integration 
time was reduced and made comparable to that for the regu- 
larized system. 

From examination of Table 4, it becomes evident that 
integration errors in the ' regularized coordinate systems 
propagate differently than do errors in the unregularized 
systems. Since a feature of regularization is the automatic 
scaling of integration step size, an increasing radius vector 
magnitude will automatically increase the step size whereas 
a decreasing radius vector magnitude will automatically 
decrease the integration step size. Thus, due to the nature 
of the Earth escape spiral trajectory, the radius vector is 
continually increasing, and it is conceivable that the step 
size will have to be reduced in order to maintain the desired 
accuracy. From examination of Table 4, it is evident that 
with only one exception, the integration step size for the 
regularized variables is always halved. The exception occurs 
for the 10 -4 to 10~ 6 error limits using the polar coordinates. 
In this case, the error is such that the 10 6 boundary is 
just crossed thereby doubling the step size. With further 
integration, the error becomes large and the step size is 
halved again. In all other instances, the lower boundaries 


11 



are never encountered. Since the lower boundaries are not 
encountered, increasing the error-bound separation limit does 
not affect the regularized systems and only penalizes the 
unregularized system by increasing the integration times. 

An alternative approach to regularization is suggested 
by the lack of encounters at the lower boundaries for the 
regularized variables. Since only the upper boundary is 
encountered, a value of n < 3/2, in the transformation 
dr = r -n dt , could be selected. This would keep the step 
size from increasing so rapidly with increasing values of 
the radius and thus eliminate the decrease in step size 
associated with an encounter with the upper boundary. Such 
a value of n would not eliminate the mathematical singulari- 
ties; however, in most normal cases the singularities are 
never encountered anyway. This concept presents an interesting 
possibility for numerical integration step size control. 

All information presented thus far has been associated 
with the characteristics of the last trajectory generated by 
an iteration process, that is the converged trajectory. It 
is of interest to know how the four different cases studied 
are affected by making certain errors in the initial assump- 
tion for boundary conditions (the Lagrange multipliers and 
terminal time) . Table 5 presents information on the number 
of iterations required and the computer time expended in 
converging from certain specified initial error percentages 
in the Lagrange multipliers. Since all possible combinations 
of the four multipliers and percentage errors represent too 
many cases to examine efficiently, all multipliers were con- 
sidered to be in error by the same percentage for each case 
studied. 
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The results presented in Table 5 indicate that the 
polar coordinates are less sensitive than the rectangular 
coordinates to errors in the initial Lagrange multipliers. 
Table 5 also indicates that regularized variables are less 
sensitive than the unregularized variables to erroneous 
initial conditions. Although the number of iterations 
required to achieve convergence is essentially the same for 
all cases, the computer time requirements are not. The 
reason that the regularized variables require less computer 
time than the unregularized variables may be seen readily by 
examining Figure 3. 

Figure 3 shows that the convergence rate of the regu- 
larized variables for initial multiplier errors of 8 percent 
is greater than the respective rate of the unregularized 
variables. The trend presented in Figure 3 is considered 
to be representative of all cases given in Table 5. Had 
Table 5 been expanded to include errors greater than ±20 
percent, the computer time savings of the regularized 
variables would probably have been more significant. Note 
that for results presented in Figure 3 and Table 5, the 
value of the terminal time was not perturbed. This, in 
general, is not realistic. If the problem is such that the 
radius vector increases with time and regularized variables 
are being used, care must be taken in the initial assumption 
for the terminal time. The sensitivity of the terminal 
pseudo time t to errors in the terminal time t in seen 
in Fig. 2. One solution involves continuously monitoring 
the terminal norm and selecting the terminal time which 
corresponds to the minimum norm for the first assumption. 
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Although for some cases the regularized and unregularized 
systems may exhibit nearly equal integration times, the inte- 
gration accuracy of each system may differ. Since a closed-form 
solution to the problem considered here does not exist, the 
error generated by the numerical integration process is 
unknown. However, there does exist a constant of motion 
which may be considered in evaluating the accuracy of the 
numerical integration procedure. This constant of motion, 
evaluated at the final time, is given by Equation 5. For 
the example cliscussed, this constant, referred to as 1+H, 
must be zero throughout the trajectory. Thus, the deviation 
of 1+H from zero is one indication of the inaccuracy of the 
numerical integration process. It should be noted, however, 
that the satisfaction of 1+H =0 is necessary but is not 
sufficient to insure numerical integration accuracy. Since 
some of the terms in the expression for 1+H contain combina- 
tions of the integrated variables, large error generation 
in two separate terms could cancel, leaving the impression 
that numerical accuracy had been achieved. 

The relative values of 1+H for converged iterations 
using the regularized and unregularized systems may be seen 
by comparing Figures 4 and 5. Figure 4 shows that the error 
in 1+H for the unregularized polar system is less than the 
error in 1+H for the rectangular system. Figure 5 indicates 
that the error in 1+H for the regularized polar system is 
larger than the error in 1+H for the regularized rectangular 
system. However, at the terminal time, the polar coordinate 
error is less than the rectangular coordinate error. Note 
also that the error in 1+H for the regularized polar system 
is quite constant during most of the integration interval; 
hence the automatic step -size adjustment associated with the 
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regularized variables tends to control the numerical error. 
Figure 4 illustrates that, for the unregularized variables, 
the error passes from a relatively large value to a relatively 
small value during the course of the trajectory. 

CONCLUSIONS 

Based on the results obtained in this study, the fol- 
lowing general conclusion can be drawn. Care in the selec- 
tion of the coordinate system used to describe an optimal 
trajectory can lead to increased accuracy and reduced 
computation time. In addition, for space vehicles subjected 
to a continuous thrust force which undergo wide variations 
in the gravitational force magnitude, significant reductions 
in computing time can be achieved by using a regularized 
form for the equations regardless of the error-bound magni- 
tude employed. In this study, reductions in computing time 
by a factor of three are obtained in some cases by using 
regularized variables. In addition, if the Hamiltonian is 
used as an indication of numerical accuracy, the trade-off 
between integration time and integration accuracy is 
apparent. It is shown that regularizing results in an 
automatic step-size change that produces relatively constant 
numerical error over the trajectory interval. These results 
indicate the importance of obtaining more definitive methods 
for selecting regularization schemes. 
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TABLE 1.- NUMERICAL INTEGRATION CHARACTERISTICS FOR ERROR BOUND SEPARATION OF 10 6 

FOR T-HE OPTIMAL LOW THRUST EARTH -ESCAPE SPIRAL 



A1 lovable 
Error 
(Absolute) 

Unregularized 

Regularized 

Rectangular 

•PolaT 

Rectangular 

Polar 

Computation time for 

10 “ 4 - io“ 10 

19 . 5 

20.6 

8.3 

7.7 

integration of state 

10 -5 - icT 11 

38.0 

21.0 

15.2 

8.1 

and perturbation 






equations (Seconds) 

10 6 - 10 

. 

71.1 

42.5 

29.4 

IS .6 

Mean computation 


■ 




time per integration 


. 0275 

. 0300 

. 0304 

.0307 

step (Seconds) 






Number of 

10" 4 - 10“ 10 

702 

685 

272 

251 

integration steps 

10~ 5 - IO -11 

1381 

702 

497 

261 


10~ 6 - 10~ 12 

2594 

1403 

971 

508 

Number of step 

io" 4 - 10 -10 

0 

1 

1 

1 

size changes 

10~ 5 - io" 11 

2 

0 

2 

2 


10“ 6 - IO -12 

3 

1 

2 

2 

Terminal error 

IO -4 - IO -10 ' 

.1375 E -10 

.4365 E -13 

.6228 E -11 

.9087 E -12 

norm 

io -5 - io -11 

.1524 E -11 

.3681 E -13 

.9438 E -09 

.8325 E -12 


io -6 - io" 12 

.2010 E -11 

.5336 E -09 

.1330 E -08 

.2150 E -11 































TABLE 2.- NUMERICAL INTEGRATION CHARACTERISTICS FOR ERROR BOUND SEPARATION OF IQ 4 

FOR THE OPTIMAL LOW THRUST EARTH ESCAPE SPIRAL 



Allowable 

Error 

(Absolute) 

Unregularized 

Regularized 

Rectangular 

Polar 

Rectangular 

Polar 

Computation time for 

10“ 4 - 10 -8 

16.4 

13.9 

8.4 

7.7 

integration of state 

10 -5 - 10 -9 

27.8 

18.2 

15.2 

8.1 

and perturbation 

a in 





equations (Seconds) 

10 b - 10 xu 

SI. 2 

31.8 

30.1 

15.7 


— 7 -11 

10 - 10 

64.0 

37.7 

34.0 

21.7 


10~ 8 - 10~ 12 

108.6 

72.4 

60.1 

32.1 

Mean computation 






time per integration 


.0276 

.0299 

.0307 

.0310 

step (Seconds) 






Number of 

10" 4 - 10~ 8 

585 

460 

272 

251 

integration steps 

10 -5 - 10~ 9 

993 

606 

497 

261 


10" 6 - 10 -1 ° 

1862 

1080 

971 

508 


10' 7 - 10 -11 

2327 

1254 

1088 

709 


10" 8 - 10 -12 

3957 

2417 

1991 

1049 




TABLE 2.- NUMERICAL INTEGRATION CHARACTERISTICS FOR ERROR BOUND SEPARATION OF 10 
FOR THE OPTIMAL LOW THRUST EARTH ESCAPE SPIRAL (Concluded) 


— — - 

Allowable 

Unregularized 

Regularized 


Error 

(Absolute) 

Rectangular 

Polar 

Rectangular 

Polar 

Number of step 

10" 4 - 10" 8 

2 ' 

2 

1 

1 

size changes 

10“ 5 - 10“ 9 

. 3 

1 

2 

2 


10~ 6 - ID -10 

4 

3 

■2 

2 


10- 7 - 10- 11 

4 

2 

3 

3 


10~ 8 - 10~ 12 

S 

3 

4 

4 

Terminal error 

10" 4 - 10“ 8 

.5603 E -10 

.1265 E -10 

.6228 E -11 

.9087 E -12 

norm 

10~ 5 - 10~ 9 

.1849 E -11 

.5304 E -13 

.9438 E -09 

.8325 E -12 


io" 6 - io" 10 

.1766 E -11 

.5328 E -09 

.1330 E -08 

.2510 E -11 


10" 7 - io' 11 

.,1413 E -11 

.5336 E -09 

.1244 E -08 

.2406 E -11 


io' 8 - io“ 12 ' 

.1378 E -11 

.6035 E -09 

.1258 E -08 

.2042 E -11 




TABLE 3.- NUMERICAL INTEGRATION CHARACTERISTICS FOR ERROR BOUND SEPARATION OF 10 

FOR THE OPTIMAL LON THRUST EARTH ESCAPE SPIRAL 


Allowable 

Error 

(Absolute) 



Computation time for 
integration of state 
and perturbation 
equations (Seconds) 


Mean computation 
time per integration 
step (Seconds) 


Number of 
integration steps 




10" 4 - 10" 6 
10" 5 - 1G~ 7 
IQ” 6 - 10' 8 

ie " 7 - ier 9 

1Q" S - 10 ~ 1 ° 
10~ 9 - 10 -11 
lO" 10 - lO ” 12 


Unregularised 

Regularized 

Rectangular 

Polar 

Rectangular 

Polar 

9.4 

7.5 

8.3 

6.1 

17.3 

10.6 

15.4 

8.1 

26.6 

15. S 

30.1 

15.7 

36.4 

26.3 

33.8 

21.7 

66.8 

43.6 

61.6 

32.6 

105.5 

60.7 

119.1 

61.2 

147.1 

102.5 

132.7 

77.8 

.0279 

.0301 

.0307 

.0307 

332 

241 

272 

193 

611 

345 

497 

261 

954 

514 

971 

SOS 

1314 

869 

1088 

709 

242'3 

1363 

1991 

1049 

3757 

2039 

3884 

2038 

5235 

3467 

4355 

2582 
















TABLE 3.- NUMERICAL INTEGRATION CHARACTERISTICS FOR ERROR BOUND SEPARATION OF 10 2 
FOR THE OPTIMAL LOW THRUST EARTH ESCAPE SPIRAL (Concluded) 



Allowable 

Error 

(Absolute) 

Unregul arized 

Regularized 

Rectangular 

. 

Polar 

Rectangular 

Polar 

Number of step 

10 -4 - 1 ( T 6 

3 

3 

1 

3 

size changes 

10~ 5 - 1 0“ 7 

4 

3 

2 

2 


10 -6 - 10 -8 

6 

4 

2 

2 


10“ 7 - 10" 9 

S 

3 

3 

3 


10 -8 - 10“ 10 

6 

S 

4 

4 


10" 9 - 10" 11 

8 

6 

4 

5 


io ” 10 - 10~ 12 

7 

5 

5 

5 

Terminal error 

10“ 4 - 10“ 6 

.2197 E -08 

. 97S0 E -13 

.6228 E -11 

.1527 E -13 

norm 

10“ 5 - 10 -7 

.1515 E -10 

.1676 E -08 

.9438 E -09 

.8325 E -12 


10" 6 - 10“ 8 

.1826 E -10 

.2231 E -09 

.1329 E -09 

.2150 E -11 


10" 7 - 10“ 9 

.2580 E -11 

.5122 E -09 

.1244 E -08 

.2406 E -11 


10" 8 - IO" 10 

.1133 E -11 

.5962 E -09 

.1258 E -08 

.2042 E -11 


10~ 9 - 10 -11 

.1624 E -11 

.6061 E -09 

.1260 E -08 

.2054 E -11 


IO -10 - IO -12 

.1560 E -10 

.6081 E -09 

.1259 E -08 

.2005 E -11 




TABLE 4. -INTEGRATION ERROR BOUNDARY ENCOUNTERS FOR VARIOUS ERROR BOUND 
SEPARATIONS FOR THE OPTIMAL LON THRUST EARTH ESCAPE SPIRAL 










TABLE 5.- INITIAL ERROR INFLUENCE ON THE CONVERGENCE CHARACTERISTICS FOR 


UNREGULARIZED AND REGULARIZED RECTANGULAR AND POLAR COORDINATES 
FOR INTEGRATION ERROR BOUNDS OF 10 _s TO 10" 9 


Initial 
Error 
In X 

Unregular ized 

Regularized 

Rectangular 

Polar 

Rectangular 

Polar 

Iterations 
Required For 
Convergence 

Computation 
Time (min) 

Iterations 
Required For 
Convergence 

Computation 
Time (min) 

Iterations 
Required For 
Convergence 

Computation 
Time (min) 

Iterations 
Required For 
Convergence 

Computation 
Time (min) 

+ 20 

6 

2.9 

5 

1.5 

6 

1.7 

5 

0.8 


5 

2.3 

5 

1.5 

6 

1.7 

5 

0.8 


S 

2.4 

4 

1.1 

5 

1.4 

4 

0.6 

+ 8 

5 

2.4 

4 

1.1 

5 

1.4 

4 

0.6 

mm 

4 



1.8 

4 

1.1 



5 

1.4 

4 

0.6 

0 

0 

0.06 

0 

0.04 

0 

0.04 

0 

0.03 

- 4 

5 

2.3 

4 

1.2 

5 

1.7 

4 

0.6 

- 8 

6 

2.9 

4 

1.2 

6 

1.7 

4 

0.6 

-12 

9 

4.7 

4 

1.2 

13 

4.2 

4 

0.6 

-16 

HKMi 

3.5 

4 

1.1 

6 

1.7 

4 

0.6 

-20 

Warn 

3.5 

4 

1.1 

6 

1.7 

5 

0.7 











































EARTH RADII 


+4 



I 1 1 l I I 

-6 -4 -2 0 +2 +4 

X, EARTH RADII 


Figure 1. - Optimal low thrust Earth escape 
spiral trajectory for T/M =0.1 
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Figure 2. - Real time vs regularized time for 
the opti mal low thrust Earth escape spiral trajectory 
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Figure 3. - Terminal error norm vs computational time 

for a SA 0 = + 8% and dty = 0 
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Figure 4. - Error in 1+H for the unregularized 
rectangular and polar coordinates for an error hound 
of 10 -5 to "SO -9 (rectangulars took 993 steps and 

polars took 606 steps) 
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Figure 5. - Error in 1 + H for the regularized 
rectangular and polar coordinates for an error bound 
of 10 -5 to 10-*(rectangulars took 497 steps 
and polars took 261 steps) 







APPENDIX A 

RECTANGULAR COORDINATES - UNREGULARIZED 



RECTANGULAR COORDINATES - UNREGULARIZED 


The equations of motion for the unregularized rectangular 
coordinates are: 


u 

= 

TA 

TJX u 

r 3 MA 

• 

V 

= 

TA 

MX v 

3 MA 

r 

• 

X 

= 

u 

• 

y 

= 

V 

M 

= 

- 6 

where 



r 


J 2 ^ 2 

Vx + y 

A 

= 

Va 2 + A 2 

U V 

U 

= 

gravitational constant 

T 

= 

thrust 

e 

— 

mass flow rate 


A-l 



The Eulex-Lagrange equations are; 


A = - w 

u u 


A 


v 


(I) 

u 


- Q1 

V 

yA u 3y(xA u + yA y )x 


0 ) 


uA v 3y(xx u + y^ v )y 


M 


TA 

„2 

M 


A- 2 



The coefficients for the perturbation equations (nonzero 
terms) are: 


9u 

3x 



3yx 2 
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* 

3u _ 3yxy 
37 x 5 


3u 

3M 


TA 

u 

M 2 X 


3u 

3A 
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3u 
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u v 

MA 3 


3v _ 3yxy 
3x r 5 


3v 

3y 



3v 

3M 


TA 
v 

M 2 A 


3v 
3 A 

u 


TA A 
u v 
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T x 2 

1 Z - 1 0 
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3 CD 

\ 
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v 3 yxy 
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2 TA 

M 3 


3 A 
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TA 

u 
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v 

M 2 A 
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The terminal 

boundary conditions in the unregularized 

rectangular 

coordinates are : 

H i * 

0.5(i 2 * y 2 ) - H. 


3 

r w u 

H_ - 

X u 

2 

u yx 


3 

r w v 

H = 

X - u 

3 

v yx 

H . = 

w y 

u' 

4 

V X 

H = 


5 

M 

cn 

ii 

i.° - ^ Cxx u + yX„) ' jt + uoi u + v«i v - X M S 
r 



The time derivatives of the terminal 


constraints are: 


h * uu + vv + ( Ux + vy) 


h = x _ 3 rUM u (ux * v yJ _ A ; r \, u t 
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px 


H 


. _ 3rv^ (u x * vy) r 3 vj T \ j r 3 a 

v ~ - _iL_ + 


UX px 2 
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v X 


V v r W ,J U 
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M 
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3H 


The nonzero elements of the mat 
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APPENDIX B 


RECTANGULAR COORDINATES - REGULARIZED 



RECTANGULAR COORDINATES - REGULARIZED 


The equations of motion for the regularized rectangular 


coordinates 

are : 


u* 


_ . 3(ux + vy) u 

Tr 3 X 

U 


2r 2 

MX 

v' 

= 

- yy + 3(ux + vy)v 
2r 2 

Tr 3 X 

V 

MX 

x> 

= 

u 


y* 

= 

V 


M' 

= 

- » r ! 


where 




r 

” 

Vx 2 + y 2 


X 

= 

Jx 2 + X 2 

T u v 


U 

= 

gravitational constant 

T 

= 

thrust 


0 

= 

mass flow rate 
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The Euler-Lagrange equations axe 
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? 
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The coefficients for the per 
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The terminal boundary conditions in the regularized 
rectangular coordinates are: 
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The time derivatives of the terminal constraints are: 
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APPENDIX C 


POLAR COORDINATES - UNREGULARIZED 



POLAR COORDINATES - UNREGULARIZED 


The equations of motion for the unregularized polar 
coordinates are: 
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where 


p = radius 
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p = gravitational constar 
T =■ thrust 
3 = mass flo\rf rate 
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The Euler-Lagrange equations are 
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The coefficients for the perturbation equations (nonzero 
terms) are: 
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The terminal boundary conditions in unregularized polar 
coordinates are: 
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The time derivatives of the terminal constraints are 
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The nonzero 
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of the matrix are; 
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APPENDIX D 

POLAR COORDINATES - REGULARIZED 



POLAR COORDINATES - REGULARIZED 


The equations of motion for the regularized polar 
coordinates are: 
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where 

p =* radius 
X - Jx 2 + X 2 

? H V 

U = gravitational constant 
T = thrust 
3 = mass flow rate 
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The Euler-Lagrange equations axe 
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he coefficients for the perturbation equations (nonzero 
terms) are: 
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The terminal boundary conditions in the regularized polar 
coordinates are: 
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The time derivatives of the terminal constraints are 
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The nonzero elements of the matrix 
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APPENDIX E 


NORMALIZED VALUES 



NORMALIZED VALUES 


In order to enhance the numerical integration accuracy, all 
numerical calculations were made in a normalized system. 

The units of normalization are given in Table E-l. The 
unit of length corresponds to one Earth radius and the unit 
of velocity to the circular velocity at one Earth radius . 

The unit of mass was chosen to be 5000 kg. The remaining 
normalization units are' such that consistent dimensional 
properties are maintained. 

Table E-2 gives the normalized values of the constants 
common to all of the coordinate systems investigated. 

Since these constants are normalized, the units are 
indicated by the general notation of L for length, 

T for time, and M for mass. 

Tables E-3 and E-4 present, respectively, the normalized 
values of the initial and terminal states for all coordinate 
systems investigated. Again, the dimensions are indicated 
by the general notation. 
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TABLE E-l. -NORMALIZATION UNITS 


Unit 

Value 

Length 

0.63781450 x 10 7 m 

Velocity 

.79053881 x 10 4 m/sec 

Time 

.80680985 x 10 3 sec 

Mass 

.5000 x 10 4 kg 

Force 

.48991644 x 10 5 Og-m)/sec 2 


TABLE E- 2 . - NORMALIZED VALUES OF CONSTANTS 


Constant 

Value 

Thrust 

Mass flow rate 
Gravitation 

0.10205822 x 10 _1 ML/T 2 
.16336057 x 10“ 5 M/T 
.10 x 10 1 L 3 /T 2 
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TABLB E-3. - NORMALIZED INITIAL CONDITIONS 


Variable 

Rectangular 

Polar 

Unregularized 

Regularized 


Regularized 

TIME (T) 

0.0 

0.0 

O 

o 

0.0 

u(L/TD 

0.0 

0.0 

0.0 

0.0 

v (L/T) 

0.97728258 

0.10470436X10 1 

0 .97728258 

0.10470436X10 1 

x(L) 

0.1047039SX10 1 

0.10470395X10 1 

0.10470395X10 1 

0.1047039SX10 1 

y (L) 

0.0 

0.0 

0.0 

0.0 

M (M) 

1.0 

1.0 

1.0 

1.0 

x u (t 2 /l) 

0.29606237X10 1 

0.2960491X10 1 

0.29608441X10 1 

0.2960I179X10 1 

X v (T 2 /L) 

-0.9792 80 73xl0 2 

-0 . 97 91 7 3 9x1 0 2 

-0.97927892xl0 2 

-0.97975524X10 2 

u u (T/L) 

-0.95538761X10 2 

-0.10234806xl0 3 

-0.95538506X10 2 

-0.10240S78xl0 3 

W v (T/L) 

0.27633966X10 1 

0.29604389X10 1 

0.27635833X10 1 

0.29607177X10 1 

x h (t/m) 

0.78700772xl0 2 

0.78697428xl0 2 

0 . 78700659x 10 2 

O.78709925xlO 2 


TABLE E-4. - NORMALIZED TERMINAL CONDITIONS 


Variable 

Rectangular 

Polar 

Unregularized 

Regularized 

Unregularized 

Regularized 

TIME(T) 

0.70l4S389xl0 2 

0. 2306334 Sxio 2 

0.7014S336xl0 2 

0.23063301xl0 2 

UCL/T) 

0.26064303 

0.64876389X10 1 

0.30879017 

0.76866S63X10 1 

v(L/T) 

-0.40823787 

-0.10162287xl0 2 

0.37315096 

0.92887282X10 1 

x(l0 

-0. 26111336X10 1 

-0.26114617X10 1 

0.8S254035X10 1 

0.85254079X10 1 

y(D 

-0.81156958X10 1 

-0.81154810X10 1 

0 . 232 S0630xl0 2 

0.23250559X10 2 

M(M) 

0.99988541 

0.99988541 

0.99988541 

0 .99988541 

1 U CT 2 /L) 

-0. S2721878xl0 2 

-0 . 5 2718636x 10 2 

-0.62460890xl0 2 

-0.62461087xl0 2 

x v Ct 2 /l) 

0 . 82 S7680Qxl0 2 

0.82578870X10 2 

-0.75479544xl0 2 

-0,75479381xl0 2 

f^UA) 

0.85237112 

0.21220771xl0 2 

-0.27830104X10 1 

-0 .69276707X10 2 

<o v (T/LI 

0.26492650X10 1 

0.65946S01xl0 2 

-0. 1864318 6xl0” 14 

-0.3S507188X10 -12 

X m (T/M) 

0 . 2 24 2 38 3 3* IQ -1 2 

0.49770030xl0 -10 

0.14723 46 6xl0 -13 

-0.16084963xl0 -12 
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